Supersymmetric Hybrid Inflation Redux 



Mansoor Ur Rehman, Qaisar Shafi, Joshua R. Wickman 

Bartol Research Institute, Department of Physics and Astronomy, 
University of Delaware, Newark, Delaware 19716, USA 
August 2009 

Abstract 

We discuss the important role played during inflation by one of the soft supersymmetry 
breaking terms in the inflationary potential of supersymmetric hybrid inflation models. 
With minimal Kahler potential, the inclusion of this term allows the prediction for the 
scalar spectral index to agree with the value = 0.963l[j:[]i^ found by WMAP5. In 
the absence of this soft term, and by taking into account only radiative and supergravity 
corrections, it is well known that > 0.985. This same soft term has previously been 
shown to play a key role in resolving the MSSM ^ problem. The tensor to scalar ratio r 
is quite small in these models, taking on values r < 10~^ in the WMAP5 2a range of n^. 

Supersymmetric (SUSY) hybrid inflation models [H El El H] , incorporating low (~ TeV) scale 
supersymmetry and minimal (A^ = 1) supergravity (SUGRA) [5], have attracted a great deal 
of attention because of their intimate connections to mainstream particle physics [6j. Inflation 
in these models is typically realized within the framework of a supersymmetric gauge theory 
based on a gauge group G, which is spontaneously broken after inflation is over to its subgroup 
H (regular inflation), or during the inflationary phase (shifted inflation [7]). There are several 
important features of supersymmetric hybrid inflation which are worth noting here. First, the 
symmetry breaking scale M is determined by the consistency of the inflationary scenario and 
turns out to be comparable to Mqut ~ 10^^ GeV. Roughly speaking, the CMB anisotropy 
ST/T is proportional to {M/MpY, which explains why M is of order Mqut [ij- Second, 
for a wide range of parameters, the inflaton field takes values that are not much larger than 
Mgut- Thus, in contrast to chaotic inflation in which the inflaton acquires trans-Planckian 
values, the supergravity corrections in supersymmetric hybrid inflation models are adequately 
suppressed [H |9] . Third, it was shown in Ref . [9] that one of the soft supersymmetry breaking 
terms generated in minimal supergravity, can play an important role during inflation. This 
term had earlier been shown to play a key role in the resolution of the MSSM /i problem |10] . 
Finally, if G is identified with [^t/(3) x SU{2) x U{1)] x U{X)b-l or SO(IO) for example, the 
observed baryon asymmetry can be generated in these models via non-thermal leptogenesis 
[TT] . The main purpose of this letter is to carefully investigate the role during inflation of the 
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soft supersymmetry breaking term. We find that its inclusion in the infiationary potential gives 
rise to a scalar spectral index which is smaller than 0.985, the previous lower bound obtained 
in these models, and therefore in better agreement with the central value of 0.963 estimated 
by the WMAP 5 yr analysis (WMAP5) [12]. Note that this result is achieved by employing 
the canonical (minimal) Kahler potential. For a discussion involving a non-minimal Kahler 
potential, see Refs. [T3| 

SUSY hybrid infiation is defined by the superpotential [21 [I] 

W = nSi^t- M"^), (1) 

where is a gauge singlet superfield and $, $ are a conjugate pair of superfields transforming 
as nontrivial representations of some gauge group G. As a simple example, G can be the 
standard model gauge group supplemented by a gauged U{1)b-l, which requires, for anomaly 
cancellations, the presence of three right handed neutrinos. W has a U{l)Fi '-R-symmetry' such 

that W e'-°'W, S e^^S*, <l><l> -> <!>$ and it can readily be seen that Eq. is the most 
general renormalizable superpotential consistent with both i?-symmetry and G-invariance. 
The SUGRA scalar potential is given by 

Vf = e'^l'^l (^Kr.^D,^WD,^W* - Snip^ \Wf^ , (2) 

with Zi being the bosonic components of the superfields Zi G {$, $, S*, • • ■ } and where we have 
defined 

D,*W* = (D^W)* and mp = Mp/VStt ~ 2.4 x 10^^ GeV is the reduced Planck mass. The 
minimal Kahler potential can be expanded as 

K = \S\^ + \^\^ + \t\\ (3) 

In the D-fiat direction |$| = |$|, and using Eqs. ([11 13]) in Eq. ([2|), the tree level global SUSY 
potential is given by 

Vf = K^M^ -\^\y + 2K^\S\^\^\^. (4) 

Assuming suitable initial conditions, the fields get trapped in the infiationary valley of local 
minima at \S\ > Sc = M and I'I'I = |$| =0, where G is unbroken. The potential is dominated 
by the constant term Vq = k^M^, thus SUSY is broken during infiation. Infiation ends when 
the infiaton drops below its critical value Sc = M and the fields roll towards the global SUSY 
minimum of the potential \S\ =0 and |$| = |$| = M. Taking into account leading order 
SUGRA corrections, as well as radiative corrections [1] and soft SUSY breaking terms, the 
potential along the infiationary trajectory (|$| = |$| =0) is of the form 

V ^ '^'M^(l+(—X- + ^Fix) + a(^) + (^y], (5) 
I \mpj 2 87r2 ^ ^ V kM J V kM / / ' ^ ^ 
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where 
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and 



a = 



2\2-A\ cos[arg5 + arg(2 - A)]. 
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Here J\f is the dimensionahty of the representation of the fields $ and $, Q the renormahzation 
scale and x = \S\ /M. In Eq. ([5]), the F(x)-term represents the radiative corrections, while the 
last two terms are the soft SUSY breaking linear and mass-squared terms, respectively. The 
quantity A — 2 represents the coefficient of the linear soft term and the magnitude of 
A is expected to be of order unity or so. The form of these soft terms comes from gravity- 
mediated SUSY breaking [10], which we employ throughout. Relative to the linear soft term, 
the soft mass term is subdominant for \a\ > 10~^. [Note that the soft linear term induces a 
VEV proportional to 1713/2 for S, which has been exploited in Ref. [10] to resolve the MSSM fi 
problem.] By virtue of the minimal Kahler potential, an \S\'^ term is not induced by SUGRA 
corrections, as this contribution is canceled by an identical term arising from the superpotential. 
As a result, the rj problem is under control in these 'minimal' models. 

In our numerical calculations, we will take m3/2 = 1 TeV and A/" = 1 for simplicity. As 
pointed out in Ref. [9], argS* can vary significantly for k ~ 10~^, yet it is always possible to 
suppress this variation via the choice of initial conditions, as explicitly shown in Ref. [L4:\. For 
K values much different from 10~^, the variation of argS* is negligibly small. Therefore, we will 
assume that arg S is constant during inflation, and will study only the speciflc cases a = and 
a = — 1. For an analysis involving a > 0, see Refs. [9l [T6] . 

For a > 0, the potential in Eq. ([5]) monotonically increases with x. In contrast, for a < 0, 
the contribution from this term can dominate in some region of parameter space for inflation 
beginning around x ~ 1. In this case, the potential develops a local maximum, giving rise 
to 'hilltop inflation' beginning near this maximum [17]. Recent treatments have shown that 
hilltop solutions exist in these models if a non- minimal Kahler potential is used [18]; here, we 
are able obtain such solutions by employing the minimal Kahler potential. It is interesting to 
note that hilltop solutions are also generated in models of non-SUSY hybrid inflation which 
include fermion-dominated radiative corrections [TU] . 

The number of e-folds after the comoving scale / has crossed the horizon is given by 



where \Si\ = x/ M is the fleld value at the comoving scale I, and Xe denotes the value of x at 
the end of inflation. If the slow-roll approximation holds, inflation ends via a waterfall induced 
at Xe = 1, although in practice the slow-roll conditions are typically violated at values of Xg 
slightly larger than unity. During inflation, the comoving scale corresponding to ko = 27t/Iq = 
0.002 Mpc~^ exits the horizon such that the number of e-folds is approximately given by 
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(9) 



where is the reheat temperature, and the subscript '0' indicates that the values are taken at 
ko. 
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After the end of inflation, the fields fall toward the SUSY vacuum and undergo damped 
oscillations about it. The coupled S", $, $ system then decays into right-handed neutrinos 
and sneutrinos, which in turn give rise to the observed baryon asymmetry via non-thermal 
leptogenesis [HI [20]. In this case, the reheat temperature is well approximated by [9] 

where minf = \/2kM is the mass of the inflaton, and we will take the mass of the heaviest light 
neutrino to be rriy^ ^ 0.05 eV. 

The amplitude of the curvature perturbation is given by 

M ( ^3/2 



where iS.-ji = 4.91 x 10~^ is the WMAP5 normalization at ko [T2]. For Af = 1, corresponding 
say to the breaking of a U(l) gauge symmetry, cosmic strings are produced upon symmetry 
breaking. The contribution of cosmic strings to the curvature perturbation goes as the string 
tension, {6T/T)cs oc Gfi ~ {M/mp)"^ [E]. As we will see, the values of M obtained in these 
models are roughly ~ 1.1 x 10^^ GeV for rig ~ 0.963, and so the contribution from cosmic 
strings is subdominant. In our calculations, we will suppress this relatively small (< a few 
percent) contribution from cosmic strings, and consider the curvature perturbation to arise 
from primordial inflation to within the uncertainty in its measured value. 
The usual slow-roll parameters may be defined as 

ml (d^Vy ml fdlV\ m% ( d^V dlV \ 

AM^ \ V ) ' ' 2M2 V ^ ; ' ^ 4M4 \ )' ^ ' 

In the slow-roll approximation (i.e. e, |?7|,'C^ ^ 1), the spectral index is given (to leading 
order) by 



(13) 



Using Eqs. (!5l)-( jTT|) . we have calculated Ug and M as a function of k for a = — 1 and a = 0, 
as shown in Fig. [TJ (In these numerical calculations, we have used the next-to-leading order 
expressions in the slow-roll approximation for added precision.) Similar to the a > case, we 
obtain two branches of solutions for a < 0, one appearing for low values and another for high 
values of M. For a = — 1 these two branches are disconnected, whereas for a = 0, they merge 
in the region where the soft SUSY mass term is important. We find that a red spectral index 
in agreement with the WMAP5 central value, Ug = 0.963, can be generated for negative values 
of a with magnitude \a\ > 10~^. For smaller |a|, a solution coinciding with a = in the rig-K 
plane is abruptly reached, and a red tilted spectrum is not exhibited for k < 10~^. 

The slow-roll parameter e is quite small in these models, and its contribution to (n^ — 1) is 
masked by the r] term. Thus we have suppressed this term in the right-hand side of Eq. (fT3!) . 
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Figure 1: Us and logiQ(M/GeV) vs. log^^o ^i "^3/2 = 1 TeV, a = and a = — 1, with a given 
by Eq. ([7]). The left panel includes the WMAP5 68% and 95% confidence level contours for the 
value of the spectral index [12j. 



The form of this equation suggests that (n^ — 1) receives a positive contribution from all terms 
except the radiative correction term, which turns out to be negative. However, the role played 
by the a-term is not readily apparent from this approximate equation; this term contributes 
most directly through the A?^ constraint (encoded via a;o), and it can be seen from Fig. [T]that 
its inclusion can lead to qualitatively different results. We now turn to a more quantitative 
analysis of the role of the a-term. 

From Eq. ([5]), we see that the SUGRA term dominates the potential if a; ^ 1. In other 
words, the limit Xq ~ 1 is a good approximation in regions where the terms in the potential 
may compete with one another. In this limit, Eq. flTT]) becomes 



A. 
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(14) 



By fixing At^ at the WMAP5 value and differentiating, it is possible to find limiting values of k 
and M with respect to one another in regions where two of the terms in Eq. (jH]) compete. As 
seen in Fig. [H each case for a will lead to a limiting for each of k and M (on separate branches 
in the a < case). 

The limiting values of k, may be obtained by differentiating Eq. f lT4|) with respect to M and 
setting dn/dM = 0. For a = 0, k is small in the limiting region, and the contribution from 
radiative corrections is negligible. The interaction of the SUGRA and SUSY mass terms then 
leads to 



Ka=0 > 



55 



m 



3/2 



mp 



9.1 X 10" 



with M ^ 2.4 X 10^^ GeV 



(15) 



In the a = —1 case, the SUSY mass term is negligible as compared to the a-term, and can be 
neglected throughout. In the limiting region of k, the SUGRA term is subdominant, and the 
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Figure 2: log;^o'^ logio \dns/dliak\ vs. Ug, for 1713/2 = 1 TeV, a = and a 
given by Eq. ([7]). 
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competition of the radiative correction and a-terms yields 



231/2^7 



37/2(AAln2)4A7^ \ mp J 



\a\m3/2 \' 



1/10 



3.4 X 10-^ with M ^ 2.0 X 10'^ GeV. (16) 



Similarly, the limiting values of M are obtained by taking dM/dn = 0. For a = 0, k is larger 
than its limit and M is at its smallest value, resulting in suppression of the SUGRA term. The 
interplay of the radiative correction and SUSY mass terms then gives 



Ma=o > mp 



n 



)3/2 "^3/2 

mp 



1/7 



.4 X 10'* GeV, with k ^ 4.7 x 10"^ (17) 



Finally, the limiting value of M appears only on the upper branch of results in the a = — 1 
case, where large M and small k values lead to SUGRA dominance over the radiative correction 
term. The interplay of the SUGRA and a-terms leads to 



-1 ^ 



\a\m3/2 



.A/67rA7j mp 



1/6 



2.6 X 10'^ GeV, with k ^ 6.0 x 10" 



(Note that while the numerical values in Eqs. ( TT6l) and ( ITSl) are given for a = — 1, we find 
numerically that the analytical expressions hold for a < — 10~^, where the soft SUSY mass 
term begins to play a significant role.) For a fixed value of k, (M) larger than its limit, there 
exist two values of M (n) which satisfy the A-ji constraint. 

The tensor to scalar ratio r is an important cosmological parameter, as it describes the 
(squared) amplitude of primordial gravitational waves and gives information on the energy 
scale of infiation. To leading order in slow-roll, this quantity is given by 
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Figure 3: logiQiTj. / GeV) vs. n and log^^ol^/GeV) vs. log]^o("^mf/GeV) for 7713/2 = 1 TeV, 
a = 0, a = —1 and a = +1 (right panel), with a and T,. given by Eqs. ([7]) and (fTOj) . respectively. 



where the first derivative of F{x) is always positive during inflation. An accurate measurement 
of r can be quite powerful in discriminating between various classes of inflation models, and 
even between qualitatively similar models (see, for instance, Ref. [21j). 

Experimental constraints on r set by WMAP have become increasingly restrictive, and 
these bounds will become even more precise with measurements to be taken with the PLANCK 
satellite. In SUSY inflation models, r is typically quite small. In the models which we currently 
consider, r takes on values r < 10~^ in the region of Ug favored by WMAP5 (see Fig. |2]). 

Precision measurements to be performed by PLANCK may also place greater constraints 
on the running of the spectral index dus/dlnk, which we find to be favored at a value of 
\dns/dlnk\ ~ 10~^. This quantity is nonzero only for next-to-leading order in slow-roll 

^ ~ 16er/ - 24e2 - 2^'. (20) 
am 

Even so. Fig. |2] shows that \dns/dlnk\ turns out to be much larger than r, particularly in the 
vicinity of the WMAP5 central value of rig. This can be understood first by noting that |?7| 3> e. 
Moreover, it can be shown that oc ^/e, which is much larger than e in the region of interest. 

An important constraint on supersymmetric inflation models arises from considering the 
reheat temperature after inflation, taking into account the gravitino problem which requires 
that Tr < 10^-10^" GeV [22]. This constraint on depends on the SUSY breaking mechanism 
and the gravitino mass For gravity mediated SUSY breaking models with unstable 

gravitinos of mass m3/2 ^ 0.1-1 TeV, Tr < 10^-10^ GeV [23], whereas Tr < 10^° GeV for stable 
gravitinos [21]. The behavior of Tr as a function of k, is shown in Fig. |3l In the WMAP5 
favored region (for a = — 1), Tr~2x 10^ GeV. For a = 0, the reheat temperature obtains a 
lower bound T,. ~ 3 x 10^ GeV. In addition, the temperature scale of reheating should always 
be lower than the mass of the inflaton, m[^f. 

Besides the upper bound placed on Tr by thermal production, there are also constraints 
arising from non-thermal gravitino production via the direct decay of the inflaton. While these 
constraints can be rather severe for SUSY hybrid inflation [25], gravitino production depends on 
the SUSY breaking sector and these models are still viable. As displayed in Fig. |3l significantly 
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lower values of mj^f can be obtained with a < 1. This extends the allowed range of parameters 
where the gravitino constraint can be evaded. With minimal Kahler potential, it can be shown 
that even the most stringent constraint can be satisfied if a is sufficiently small. 

To summarize, we have investigated the important role played during hybrid inflation by one 
of the soft supersymmetry breaking terms generated in minimal supergravity. In the simplest 
example a scalar spectral index value of Ug = 0.963, preferred by WMAP5, is realized if the 
overall sign of this term is negative. This case also permits one to estimate the infiaton mass to 
be of order 10^ — 10^^ GeV. It is amusing that this same soft term, linear in the infiaton field, 
also allows one to resolve the MSSM /i problem. Thus we conclude that supersymmetric hybrid 
inflation models with minimal Kahler potential are both realistic and in good agreement with 
the current observations. An important prediction of these models has to do with the tensor 
to scalar ratio r. We find that r < 10"^ for 0.928 < Us < 1.008. 
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